The role of non-coding RNA in the development of pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is a progressive disease of the small pulmonary arteries, characterised by pulmonary vascular remodelling due to excessive proliferation and resistance to apoptosis of pulmonary artery endothelial cells (PAECs) and pulmonary artery smooth muscle cells (PASMCs). The increased pulmonary vascular resistance and elevated pulmonary artery pressures result in right heart failure and premature death. Germline mutations of the bone morphogenetic protein receptor-2 (bmpr2) gene, a receptor of the transforming growth factor beta (TGF-β) superfamily, account for approximately 75%-80% of the cases of heritable form of PAH (HPAH) and 20% of sporadic cases or idiopathic PAH (IPAH). IPAH patients without known bmpr2 mutations show reduced expression of BMPR2. However only ~ 20% of bmpr2-mutation carriers will develop the disease, due to an incomplete penetrance, thus the need for a ‘second hit’ including other genetic and/or environmental factors is accepted. Diagnosis of PAH occurs most frequently when patients have reached an advanced stage of disease. Although modern PAH therapies can markedly improve a patient’s symptoms and slow the rate of clinical deterioration, the mortality rate from PAH remains unacceptably high. Therefore, the development of novel therapeutic approaches is required for the treatment of this multifaceted disease. Noncoding RNAs (ncRNAs) include microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). MiRNAs are ~ 22 nucleotide long and act as negative regulators of gene ex-pression via degradation or translational inhibition of their target mRNAs. Previous studies showed extensive evidence for the role of miRNAs in the development of PAH. LncRNAs are transcribed RNA molecules greater than 200 nucleotides in length. Similar to classical mRNA, lncRNAs are translated by RNA polymerase II and are generally alternatively spliced and polyadenylated. LncRNAs are highly versatile and function to regulate gene expression by diverse mechanisms. Unlike miRNAs, which exhibit well-defined actions in negatively regulating gene expression via the 3’-UTR of mRNAs, lncRNAs play more diverse and unpredictable regulatory roles. Although a number of lncRNAs have been intensively investigated in the cancer field, studies of the role of lncRNAs in vascular diseases such as PAH are still at a very early stage. The aim of this study was to investigate the involvement of specific ncRNAs in the development of PAH using experimental animal models and cell culture. The first ncRNA we focused on was miR-143, which is up-regulated in the lung and right ventricle tissues of various animal models of PH, as well as in the lungs and PASMCs of PAH patients. We show that genetic ablation of miR-143 is protective against the development of chronic hypoxia induced PH in mice, assessed via measurement of right ventricular systolic pressure (RVSP), right ventricular hypertrophy (RVH) and pulmonary vascular remodelling. We further report that knockdown of miR-143-3p in WT mice via anti-miR-143-3p administration prior to exposure of mice to chronic hypoxia significantly decreases certain indices of PH (RVSP) although no significant changes in RVH and pulmo-nary vascular remodelling were observed. However, a reversal study using antimiR-143-3p treatment to modulate miR-143-3p demonstrated a protective effect on RVSP, RVH, and muscularisation of pulmonary arteries in the mouse chronic hypoxia induced PH model. In vitro experiments showed that miR-143-3p overexpression promotes PASMC migration and inhibits PASMC apoptosis, while knockdown miR-143-3p elicits the opposite effect, with no effects observed on cellular proliferation. Interestingly, miR-143-3p-enriched exosomes derived from PASMCs mediated cell-to-cell communication between PASMCs and PAECs, contributing to the pro-migratory and pro-angiogenic phenotype of PAECs that underlies the pathogenesis of PAH. Previous work has shown that miR-145-5p expression is upregulated in the chronic hypoxia induced mouse model of PH, as well as in PAH patients. Genetic ablation and pharmacological inhibition (subcutaneous injection) of miR-145-5p exert a protective against the de-velopment of PAH. In order to explore the potential for alternative, more lung targeted delivery strategies, miR-145-5p expression was inhibited in WT mice using intranasal-delivered antimiR-145-5p both prior to and post exposure to chronic hypoxia. The decreased expression of miR-145-5p in lung showed no beneficial effect on the development of PH compared with control antimiRNA treated mice exposed to chronic hypoxia. Thus, miR-143-3p modulated both cellular and exosome-mediated responses in pulmonary vascular cells, while the inhibition of miR-143-3p prevented the development of experimental pulmonary hypertension. We focused on two lncRNAs in this project: Myocardin-induced Smooth Muscle Long noncoding RNA, Inducer of Differentiation (MYOSLID) and non-annotated Myolnc16, which were identified from RNA sequencing studies in human coronary artery smooth muscle cells (HCASMCs) that overexpress myocardin. MYOSLID was significantly in-creased in PASMCs from patients with IPAH compared to healthy controls and increased in circulating endothelial progenitor cells (EPCs) from bmpr2 mutant PAH patients. Exposure of PASMCs to hypoxia in vitro led to a significant upregulation in MYOSLID expres-sion. MYOSLID expression was also induced by treatment of PASMC with BMP4, TGF-β and PDGF, which are known to be triggers of PAH in vitro. Small interfering RNA (siR-NA)-mediated knockdown MYOSLID inhibited migration and induced cell apoptosis without affecting cell proliferation and upregulated several genes in the BMP pathway in-cluding bmpr1α, bmpr2, id1, and id3. Modulation of MYOSLID also affected expression of BMPR2 at the protein level. In addition, MYOSLID knockdown affected the BMP-Smad and BMP-non-Smad signalling pathways in PASMCs assessed by phosphorylation of Smad1/5/9 and ERK1/2, respectively. In PAECs, MYOSLID expression was also induced by hypoxia exposure, VEGF and FGF2 treatment. In addition, MYOSLID knockdown sig-nificantly decreased the proliferation of PAECs. Thus, MYOSLID may be a novel modulator in pulmonary vascular cell functions, likely through the BMP-Smad and –non-Smad pathways. Treatment of PASMCs with inflammatory cytokines (IL-1 and TNF-α) significantly in-duced the expression of Myolnc16 at a very early time point. Knockdown of Myolnc16 in vitro decreased the expression of il-6, and upregulated the expression of il-1 and il-8 in PASMCs. Moreover, the expression levels of chemokines (cxcl1, cxcl6 and cxcl8) were sig-nificantly decreased with Myolnc16 knockdown. In addition, Myolnc16 knockdown decreased the MAP kinase signalling pathway assessed by phosphorylation of ERK1/2 and p38 MAPK and inhibited cell migration and proliferation in PASMCs. Thus, Myolnc16 may a novel modulator of PASMCs functions through anti-inflammatory signalling pathways. In summary, in this thesis we have demonstrated how miR-143-3p plays a protective role in the development of PH both in vivo animal models and patients, as well as in vitro cell cul-ture. Moreover, we have showed the role of two novel lncRNAs in pulmonary vascular cells. These ncRNAs represent potential novel therapeutic targets for the treatment of PAH with further work addressing to investigate the target genes, and the pathways modulated by these ncRNAs during the development of PAH.

Investigation of the role of long non-coding RNAs in oncogene induced cellular senescence

Cellular senescence is a stable arrest of cell proliferation induced by several factors such as activated oncogenes, oxidative stress and shortening of telomeres. Senescence acts as a tumour suppression mechanism to halt the progression of cancer. However, senescence may also impact negatively upon tissue regeneration, thus contributing to the effects of ageing. The eukaryotic genome is controlled by various modes of transcriptional and translational regulation. Focus has therefore centred on the role of long non- coding RNAs (lncRNAs) in regulating the genome. Accordingly, understanding how lncRNAs function to regulate the senescent genome is integral to improving our knowledge and understanding of tumour suppression and ageing. Within this study, I set out to investigate the expression of lncRNAs’ expression within models of senescence. Through a custom expression array, I have shown that expression of multiple different lncRNAs is up-regulated and down regulated in IMR90 replicative senescent fibroblasts and oncogene-induced senescent melanocytes. LncRNA expression was determined to be specific to stable senescence-associated cell arrest and predominantly within the nucleus of senescent cells. In order to examine the function of lncRNA expression in senescence, I selected lncRNA transcript ENST0000430998 (lncRNA_98) to focus my investigations upon. LncRNA_98 was robustly upregulated within multiple models of senescence and efficiently depleted using anti-sense oligonucleotide technology. Characterisation and unbiased RNA-sequencing of lncRNA_98 deficient senescent cells highlighted a list of genes that are regulated by lncRNA_98 expression in senescent cells and may regulate aspects of the senescence program. Specifically, the formation of SAHF was impeded upon depletion of lncRNA_98 expression and levels of total pRB protein expression severely decreased. Validation and recapitulation of consequences of pRB depletion was confirmed through lncRNA_98 knock-out cells generated using CRISPR technology. Surprisingly, inhibition of ATM kinase functions permitted the restoration of pRB protein levels within lncRNA_98 deficient cells. I propose that lncRNA_98 antagonizes the ability of ATM kinase to downregulate pRB expression at a post-transcriptional level, thereby potentiating senescence. Furthermore, lncRNA expression was detected within fibroblasts of old individuals and visualised within senescent melanocytes in human benign nevi, a barrier to melanoma progression. Conversely, mining of 337 TCGA primary melanoma data sets highlighted that the lncRNA_98 gene and its expression was lost from a significant proportion of melanoma samples, consistent with lncRNA_98 having a tumour suppressor functions. The data presented in this study illustrates that lncRNA_98 expression has a regulatory role over pRB expression in senescence and may regulate aspects of tumourigenesis and ageing.

The metabolic and epigenetic effects of the isocitrate dehydrogenase-1 mutation in glioma

Cancer cells have been noted to have an altered metabolic phenotype for over ninety years. In the presence of oxygen, differentiated cells predominately utilise the tricarboxylic acid (TCA) cycle and oxidative phosphorylation to efficiently produce energy and the metabolites necessary for protein and lipid synthesis. However, in hypoxia, this process is altered and cells switch to a higher rate of glycolysis and lactate production to maintain their energy and metabolic needs. In cancer cells, glycolysis is maintained at a high rate, even in the presence of oxygen; a term described as “aerobic glycolysis”. Tumour cells are rapidly dividing and have a much greater need for anabolism compared to normal differentiated cells. Rapid glucose metabolism enables faster ATP production as well as a greater redistribution of carbons to nucleotide, protein, and fatty acid synthesis, thus maximising cell growth. Recently, other metabolic changes, driven by mutations in genes related to the TCA cycle, indicate an alternative role for metabolism in cancer, the “oncometabolite”. This is where a particular metabolite builds up within the cell and contributes to the tumorigenic process. One of these genes is isocitrate dehydrogenase (IDH) IDH is an enzyme that forms part of the tricarboxylic acid (TCA) cycle and converts isocitrate to α-ketoglutarate (α-KG). It exists in three isoforms; IDH1, IDH2 and IDH3 with the former present in the cytoplasm and the latter two in the mitochondria. Point mutations have been identified in the IDH1 and IDH2 genes in glioma which result in a gain of function by converting α-KG to 2-hydroxyglutarate (2HG), an oncometabolite. 2HG acts as a competitive inhibitor of the α-KG dependent dioxygenases, a superfamily of enzymes that are involved in numerous cellular processes such as DNA and histone demethylation. It was hypothesised that the IDH1 mutation would result in other metabolic changes in the cell other than 2HG production, and could potentially identify pathways which could be targeted for therapeutic treatment. In addition, 2HG can act as a potential competitive inhibitor of α-KG dependent dioxygenases, so it was hypothesised that there would be an effect on histone methylation. This may alter gene expression and provide a mechanism for tumourogenesis and potentially identify further therapeutic targets. Metabolic analysis of clinical tumour samples identified changes associated with the IDH1 mutation, which included a reduction in α-KG and an increase in GABA, in addition to the increase in 2HG. This was replicated in several cell models, where 13C labelled metabolomics was also used to identify a possible increase in metabolic flux from glutamate to GABA, as well as from α-KG to 2HG. This may provide a mechanism whereby the cell can bypass the IDH1 mutation as GABA can be metabolised to succinate in the mitochondria by GABA transaminase via the GABA shunt. JMJ histone demethylases are a subset of the α-KG dependent dioxygenases, and are involved in removing methyl groups from histone tails. Changes in histone methylation are associated with changes in gene expression depending on the site and extent of chemical modification. To identify whether the increase in 2HG and fall in α-KG was associated with inhibition of histone demethylases a histone methylation screen was used. The IDH1 mutation was associated with an increase in methylation of H3K4, which is associated with gene activation. ChiP and RNA sequencing identified an increase in H3K4me3 at the transcription start site of the GABRB3 subunit, resulting in an increase in gene expression. The GABRB3 subunit forms part of the GABA-A receptor, a chloride channel, which on activation can reduce cell proliferation. The IDH1 mutation was associated with an increase in GABA and GABRB3 subunit of the GABA-A receptor. This raises the possibility of GABA transaminase as a potential therapeutic target. Inhibition of this enzyme could reduce GABA metabolism, potentially reducing any beneficial effect of the GABA shunt in IDH1 mutant tumours, and increasing activation of the GABA-A receptor by increasing the concentration of GABA in the brain. This in turn may reduce cell proliferation, and could be achieved by using Vigabatrin, a GABA transaminase inhibitor licensed for use in epilepsy.

Macrophages and CD4 T-cells in rheumatoid arthritis and their modulation by JAK inhibitors

Background: Rheumatoid arthritis (RA) is a chronic inflammatory arthritis that causes significant morbidity and mortality and has no cure. Although early treatment strategies and biologic therapies such as TNFα blocking antibodies have revolutionised treatment, there still remains considerable unmet need. JAK kinase inhibitors, which target multiple inflammatory cytokines, have shown efficacy in treating RA although their exact mechanism of action remains to be determined. Stratified medicine promises to deliver the right drug to the right patient at the right time by using predictive ‘omic biomarkers discovered using bioinformatic and “Big Data” techniques. Therefore, knowledge across the realms of clinical rheumatology, applied immunology, bioinformatics and data science is required to realise this goal. Aim: To use bioinformatic tools to analyse the transcriptome of CD14 macrophages derived from patients with inflammatory arthritis and define a JAK/STAT signature. Thereafter to investigate the role of JAK inhibition on inflammatory cytokine production in a macrophage cell contact activation assay. Finally, to investigate JAK inhibition, following RA synovial fluid stimulation of monocytes. Methods and Results: Using bioinformatic software such as limma from the Bioconductor repository, I determined that there was a JAK/STAT signature in synovial CD14 macrophages from patients with RA and this differed from psoriatic arthritis samples. JAK inhibition using a JAK1/3 inhibitor tofacitinib reduced TNFα production when macrophages were cell contact activated by cytokine stimulated CD4 T-cells. Other pro-inflammatory cytokines such as IL-6 and chemokines such as IP-10 were also reduced. RA synovial fluid failed to stimulate monocytes to phosphorylate STAT1, 3 or 6 but CD4 T-cells activated STAT3 with this stimulus. RNA sequencing of synovial fluid stimulated CD4 T-cells showed an upregulation of SOCS3, BCL6 and SBNO2, a gene associated with RA but with unknown function and tofacitinib reversed this. Conclusion: These studies demonstrate that tofacitinib is effective at reducing inflammatory mediator production in a macrophage cell contact assay and also affects soluble factor mediated stimulation of CD4 T-cells. This suggests that the effectiveness of JAK inhibition is due to inhibition of multiple cytokine pathways such as IL-6, IL-15 and interferon. RNA sequencing is a useful tool to identify non-coding RNA transcripts that are associated with synovial fluid stimulation and JAK inhibition but these require further validation. SBNO2, a gene that is associated with RA, may be biomarker of tofacitinib treatment but requires further investigation and validation in wider disease cohorts.